extern int tred1_(integer *nm, integer *n, doublereal *a, doublereal *d__, doublereal *e, doublereal *e2);
